Theoretical Priors On Modified Growth Parametrisations 
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Next generation surveys will observe the large-scale structure of the Universe with unprecedented 
accuracy. This will enable us to test the relationships between matter over-densities, the curvature 
perturbation and the Newtonian potential. Any large-distance modification of gravity or exotic 
nature of dark energy modifies these relationships as compared to those predicted in the standard 
smooth dark energy model based on General Relativity. In linear theory of structure growth such 
modifications are often parameterised by virtue of two functions of space and time that enter the 
relation of the curvature perturbation to, first, the matter over- density, and second, the Newtonian 
potential. We investigate the predictions for these functions in Brans-Dicke theory, clustering dark 
energy models and interacting dark energy models. We find that each theory has a distinct path in 
the parameter space of modified growth. Understanding these theoretical priors on the parameter- 
isations of modified growth is essential to reveal the nature of cosmic acceleration with the help of 
upcoming observations of structure formation. 
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I. INTRODUCTION 



Our understanding of the laws of physics we gained on 
Earth has already been challenged by cosmological obser- 
vations of the formation of galaxies and the accelerated 
expansion of the Universe 0, H| . In order to establish a 
cosmological standard model without inconsistencies be- 
tween observations and known physics, it is unavoidable 
either to assume the presence of as yet undetected dark 
materials filling 95% of the energy of the Universe, or 
to modify some of the fundamental laws of physics such 
as Einstein's General Relativity. The information on the 
underlying physics of the Universe is provided by look- 
ing at the formation of the large scale structure (LSS) 
of matter in the Universe. While current observations 
have revealed the breakdown of our knowledge of physics 
on cosmological scales, the future observation of the LSS 
will provide a clue to reveal the nature of the late-time 
cosmic acceleration; which part of our physical picture 
should be modified - matter or gravity, and how it is 
modified [3h32}. 

The current standard cosmological model was estab- 
lished based upon unknown dark materials and the va- 
lidity of General Relativity (GR) on cosmological scales. 
We call this standard model smooth dark energy (sDE). It 
consists of the following energy components: 5% baryons 
as described by the standard model of particle physics 
and 95% of unobserved dark particles - 20% of collision- 
less dark matter, fitting galaxy formation and the peaks 
in the spectrum of cosmic microwave background (CMB) 
anisotropics, and 75% dark energy explaining the cosmic 
acceleration, described by a perfect fluid that is homo- 
geneous (smooth) on sub-horizon scales. It is difficult to 
accept dark energy as a standard building block of energy 
components of the Universe as it is not yet supported 



neither phenomenologically nor theoretically. Gravity 
on cosmological scales can be modified to accelerate the 
cosmic expansion without dark energy. The degeneracy 
between GR and modified gravity (MG) models in the 
background expansion history is broken by the growth 
history of structure formation. Any departure from the 
growth history in sDE is probed by two functions ex- 
pressing new degrees of freedom: the mass screening, Q, 
changes the relationship between the curvature pertur- 
bation and the matter over-densities, while the effective 
anisotropic stress, 77, alters the relationship between the 
curvature perturbation and the Newtonian potential. 

Here we study whether the possible detection of non- 
trivial Q and 77 can be considered as a probe of the break- 
down of the sDE model at cosmological scales. We show 
that non-trivial Q and r) functions are also induced in 
other GR models via clustering of dark energy or inter- 
actions between dark energy and dark matter. GR will 
be confirmed (or excluded) by looking at the relation be- 
tween Q and 77, not merely by a departure of Q and/or 
77 from unity. We find that both MG models [HI, [34| 
and exotic GR models (35| provide unique distinguish- 
able trajectories that are well presented in the plane of 
two transformed parameters: S = Q(l + l/r/)/2 param- 
eterises the relation between the lensing potential and 
matter over-densities and (i = Q/r) represents the re- 
lation between the Newtonian potential and the matter 
over-densities. We examine these parameters in various 
models and show that each theory has a distinct path 
on the plane of those parameters. The understanding of 
these theoretical priors is essential to reveal the nature of 
the cosmic acceleration from future structure formation 
observations. 

This paper is organised as follows. In section II we 
introduce the smooth dark energy (sDE) model and the 
modified growth parameters Q and rj. Then we investi- 
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gate these functions for different cases: In section III the 
Brans-Dickc (BD) theory is used to characterise modified 
gravity models on sub-horizon scales. In section IV we 
consider clustering dark energy (cDE) models and in sec- 
tion V interacting dark energy (IDE) models are studied. 
In section VI we show that these models have distinct 
paths in the plane of the modified growth parameters fi 
and £. We conclude in section VII. 



II. PARAMETERISATION 

Any deviation from our standard cosmological model 
sDE will lead to modifications of the background evo- 
lution and structure growth. Let us fix the background 
using the sDE reference model and, subsequently, param- 
eterise the modified growth in terms of the functions Q 
and T] such that these become trivial (Q = r) = 1) in sDE. 

We assume a flat Friedmann- Robertson- Walker (FRW) 
background metric. In the sDE model the Friedmann 
equation defines the evolution of the Hubble parameter, 
H , in terms of the energy content of the Universe 
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(1) 



We introduced the label "s" to mark the quantities given 
in the sDE model. In any other model we adjust the 
parameters such that H evolves in the same way as in 
the sDE reference model. 

Next, we consider linear perturbations around the flat 
FRW background described in the Newtonian gauge: 

ds 2 = -(1 + 2^)dt 2 + a 2 (l + 2$)d ij dx i dx j . (2) 

We will use a dot to denote the derivative with respect 
to cosmic time, t. In GR the Newtonian potential "J and 
the curvature perturbation $ are related to the matter 
content via the linearised Einstein equations. In sDE 
these yield two constraint equations. Firstly, because 
there is no anisotropic stress at late times the two metric 
potentials coincide, ty s = — $ s , and secondly, the Pois- 
son equation relates 3> to the comoving matter density 
perturbation A m : 



fc 2 $ s = 4irG 



A? 



(3) 



For a general fluid (i) the comoving density perturbation 
is defined as 



fh(k_ 
Pi k 2 ' 



(4) 



in terms of the density contrast Si and the velocity diver- 
gence 9i in the Newtonian gauge. The comoving density 
perturbation A m is useful as it is the over-density mea- 
sured by observers comoving with the matter-flow. 

In any other GR dark energy or MG model the Poisson 
equation is modified. This departure from sDE in the 



Poisson equation is parameteriscd by the function Q(a, k) 
such that 



fc 2 $ = 4TrG N a 2 Q(a, k)p s m A r 



(5) 



holds with A m being the comoving matter perturbation 
as measured in the modified model (with the same back- 
ground evolution as sDE) and p s m that scales in the usual 
way with time (tx a -3 ). Because and A m are gauge- 
invariant variables this is a gauge-invariant definition of 
Q. In MG models, Q represents the mass screening ef- 
fect due to a local modification of gravity, i.e. an effective 
modification of Newton's constant. In exotic dark energy 
models within GR it measures the additional clustering 
due to a modification of the background evolution of the 
energy densities and different evolution of the perturba- 
tions. 

In order to fully describe the modified growth of per- 
turbations, the other constraint equation of GR needs to 
be parameterised as well. The relation between $ and 'J 
is altered such that 



(0) 



is not necessarily trivial. As compared to sDE in GR, 
where <j>/^ = —I, a non-trivial r\ effectively adds an 
anisotropic stress component at linear level. 

Measuring $ and therefore Q directly from the galaxy 
distribution is difficult due to galaxy bias. Therefore, we 
will also present our results in terms of two re-defined 
functions that are more directly connected to unbiased 
observables. First, we notice that weak lensing and the 
late ISW effect are sensitive to ($ — *)/2 and therefore 
it is convenient to use the combination 



Q 

S= - 
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I 

1 + - 

V 



(7) 



while peculiar velocities are responding to ^ and, thus, 
they determine 



/' 



V 



(8) 



Note that the growth of matter perturbations is also de- 
termined by the Newtonian potential, and hence by p. 
Therefore, these transformed functions are more closely 
related to observations. 

In the following sections we will study the evolution of 
matter fluctuations inside the Hubble horizon. Then we 
can use the quasi-static approximation and neglect the 
time derivatives of the perturbed quantities compared 
with the spatial derivatives. Furthermore, the comoving 
density perturbation and the Newtonian gauge density 
contrast are approximately equal, A^ ~ Si, which means 
we can neglect the contribution from the velocity po- 
tentials to the Poisson equation. However, this is not 
a perfect approximation when considering dark energy 
perturbations as we will discuss in Sec IIVI 
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III. QUASI-STATIC PERTURBATIONS IN 
MODIFIED GRAVITY MODELS 

As mentioned in the introduction, the large distance 
modification of gravity, which is necessary to explain the 
late-time acceleration, generally modifies gravity even on 
sub-horizon scales due to the introduction of a new scalar 
degree of freedom. This modification of gravity due to 
the scalar mode can be described by Brans-Dicke (BD) 
gravity [23|, HH, l36fl . In the quasi-static approximation 
the perturbed modified Einstein equations give 



$ + * 

',2 



k 2 

<±i:G N a 2 p m 5 m + —ip, 



(3 + 2cj B D)fc f = -8TrG N a p m S m , 



(9) 
(10) 
(11) 



where p m is the background dark matter energy density, 
S m is dark matter density contrast and if is the BD scalar 
field. In the quasi-static approximation the BD param- 
eter cjbd and Newton's constant Gjv can be any func- 
tions of time. In the following we assume that Gn can 
be approximated by the one measured in Cavendish-likc 
experiments and the it is constant. 

In general, modified gravity models that explain the 
late-time acceleration predict wbd ~ 0(1) on sub- horizon 
scales today. This would contradict the solar system con- 
straints which require wbd > 40000. However, this con- 
straint can be applied only if the BD scalar has no po- 
tential and no self-interactions. Thus, in order to avoid 
this constraint, the BD scalar should acquire some inter- 
action terms on small scales. Generally, we expect that 
the BD scalar field equation is given by [36[ 



(3 + 2u;bd)&V = —8TrGNa 2 pmS m — I(^>) 



(12) 



in a Fourier space. Here the interaction term I can be 
expanded as 



1 f d 3 k l( l 3 k 2 
(2tt) 3 

xM 2 (k 1 ,k 2 )<p(k 1 )ip(k 2 ) - 



I(tp) = M 2 <p(k) + 



S D (k - ku) x 



(13) 



where kij = ki + kj . Non-linear terms become important 
when considering non-linear clustering of dark matter 
and they are crucial to recover GR on non-linear scales. 
However, in this work we are interested in linear pertur- 
bations thus we only take into account the linear term. 

From Eqs. (j9"1)- (fT3"]) wc can directly derive the functions 
Q(a, k) and r](a, k): 



Q(k, a) 
T]{k,a) 



2(l + ^ BD ) + M 2 a 2 /fc 2 
3 + 2ujbd + M 2 a 2 /k 2 
2(1 + uj bd ) + M 2 a 2 jk 2 
2{2 + uj B D) + M 2 a 2 /k 2 



(14) 
(15) 



There are two ways to recover GR. One is to take | wbd | ^ 
1 and the other one is to consider large scales Ma/k ^> 1. 



In both cases Q and 77 approach unity. For M — 0, there 
is a unique trajectory in the Q — rj plane parameterised 
by wbd- 

As discussed in the previous section, the parameters 
that are more closely related to observations are £ and 
p; see the definitions (|7|) and (|8|). We find that these 
parameters arc given by 



S(/c, a) = 1, p(k, a) 



2(2 + w BD ) + M 2 a 2 /fc 2 
3 + 2cj B D + M 2 a 2 /fc 2 



(16) 



Interestingly, S(fc, a) = 1 is the same as in GR. This is 
because the scalar field couples to the trace of the energy- 
momentum tensor and it does not couple to the photon. 
Since £ parameterises the lensing potential, which in turn 
determines the geodesies of photons, it is not affected by 
the scalar field and, thus, £ = 1 as in GR. 

On the other hand, p measures the strength of Newton 
gravity. In order to understand the result, it is instructive 
to find a solution for the scalar filed <p with a localised 
matter source Sp m in the Mikowski background. Let us 
first consider the case with M ~ 0. Assuming spherical 
symmetry, the solution for (p(r) far away from the source 
in real space is given by 



tp(r) 



1 



3 + 2wbd r 



(17) 



where r g = 87tGjv / dr'r' 2 5p m is the gravitational length 
of the source. If 3 + 2u>bb > 0, the scalar field mediates 
an attractive force. From Eqs. © and (flU)) we obtain an 
equation for the Newtonian potential 



k 2 

k 2 ^ = -AirG N a 2 p m 5 m + —<p 



(18) 



Then we can see that the scalar field adds an additional 
attractive force to the Newtonian potential if 3 + 2wbd > 
and, consequently, p, becomes larger than one. On the 
other hand, for 3 + 2wbd < the scalar field mediates 
a repulsive force. This is related to the fact that for 
3 + 2cjbd < the scalar field becomes a ghost that has 
a kinetic term with the wrong sign. In this case, the 
scalar field suppresses the Newtonian force and p be- 
comes smaller than one. 

Next, let us discuss the effect of the mass. For 3 + 
2wbd > the solution for the scalar field in case of a 
spherically symmetric source becomes 



<p(r) 



1 



3 + 2o; B d r 



exp 



r 

tm 



where 



tm 



2^bd 



M 2 



(19) 



(20) 



is the Compton wavelength of the scalar field. On large 
scales, r > rj;, the solution exponentially decays and 
the scalar field becomes negligible. In this case there is 
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no scalar force and GR is recovered, so p, = 1. On the 
other hand, for r < rjv/, the effect of the mass can be 
neglected and we get the same result as for the masslcss 
case. Thus, for M^0,(j becomes scale-dependent. 

This BD theory contains two well known modified 
gravity models as special cases: f(R) gravity and the 
DGP braneworld model. Firstly, in f(R) gravity we have 
wed = [13]. In this case p, = 4/3 for Ma/k < 1 and 
p = 1 for Ma/k >■ 1. In f(R), though, Gjv can be time- 
dependent, which effectively modifies £ and p. However, 
its dependence on time is strongly constrained by obser- 
vations and we can safely neglect such an effect. Sec- 
ondly, in the DGP model we have w BD = (3/2) (/3(t) - I) 
with (3(t) = 1 - 2Hr c {l + H/3H 2 ), where H is the time 
derivative of the Hubble parameter and r c ~ Hq 1 is the 
cross-over scale, which is a parameter in this model [33j . 
Here M = for quasi-static perturbations. We should 
note that 3 + 2u;bd < and p < 1 in DGP indicating a 
theoretical pathology of the model [13] ■ 



IV. CLUSTERING DARK ENERGY 

The dark energy can in principle support long-lived 
fluctuations |3a - l4q ]. In these clustering dark energy 
(cDE) models the baryons will fall into the potential wells 
created by both the dark matter and energy. Here we 
consider a general dark energy fluid with a time- varying 
barotropic equation of state parameter, w de = Pde/ Pde-, 
and some scalar anisotropic stress o~de- However, in this 
section we restrict the dark energy to be minimally cou- 
pled to gravity and (dark and baryonic) matter, which 
means that here pi = p\ for i = m, de; and we can treat 
the baryons and the dark matter as a single pressureless 
perfect fluid. 

Given the different physical behaviour of matter and 
dark energy the total density fluctuation will not be a 
simple linear function of 8 m but acquires a contribution 
from the dark energy The Poisson equation in cDE mod- 
els is then given by 



k 2 $ = 4irG N a 2 (p m 5 m + p de 8 de ) 



(21) 



The dark energy anisotropic stress, o~ de , enters the second 
Einstein constraint equation as 



fc 2 ($ + = -\2irG N a 2 (l + w de )p de a de 



(22) 



Because the evolution of the background energy densities 
in cDE is not modified with respect to the sDE reference, 
we have p s m = p m . So the definition of Q in Eq. ((SJ) 
directly implies 



Q(a,k) = 1 + 



Pde$, 



A 1 
Pm^m 



(23) 



for cDE models. By combining the two constraints we 
find that the rj parameter here is 



PmSm + Pde$de 



PmSm + PdeSde + 3(1 + W de )p de a de 



(24) 



As we want to get a feeling for the time- and scale- 
dependence of Q and rj we take a look at the evolution 
equations of S m and Sde without specifying the dark en- 
ergy pressure and anisotropic stress perturbations. Later 
we will consider a simple model for the pressure pertur- 
bation and a specific model for the anisotropic stress of 
dark energy motivated by the possibility of having scale- 
invariant growth of dark energy perturbations. 

From the energy-momentum conservation of the mat- 
ter and the dark energy we can derive second-order dif- 
ferential equations for S m and 8 de in the quasi-static ap- 
proximation where time derivatives of the gravitational 
potentials are neglected. For the matter component we 
find the well-known growth equation 



8 m + 



2HS m = — > 



(25) 



with the source term as found from the two Einstein con- 
straint equations above: 

k 2 3 

= T,H 2 {^mS m + VL de 5 de + 3(1 + W de )VL de a de } ■ 

(26) 

Here fii = f2j(i) are the fractional energy densities as 
functions of time. For the general barotropic dark energy 
fluid perturbations we find the evolution equation 



S de + (2 - 6w de )H5 de + 3H 



[5P dl 



\ Pde 



3Hw de S de + 3 



-(1 + W de )- 



(2 - 3w de )H 2 + H 



w de S d 



de 



Pde 



de 



(Jde + * 



(27) 



(1 + w de )p de 

A priori there are no evolution equations for 5P de and 
o~ de available. To close the system we need some relation 
between the pressure and anisotropic stress perturbations 
and the density perturbation. 

Formally, cDE can predict any Q(a,k) and rj{a,k). 
From Eqs. (|23|) and (|24|) . we can derive 5 de and a de in 
terms of Q and 77. Then substituting them into Eq. (|27|). 
we can find a solution for 8 Pde- However, in physical 
models, 8P de and a de should be determined by micro- 
physics of the dark energy and it is very unlikely that 
we have enough degrees of freedom to fine-tune these 
quantities to get any desirable Q and rj. Of course it is a 
formidable task to find a theoretical prior for dark energy 
clustering. In the following, we consider a toy model for 
8P de and a de to get some idea. 

The pressure perturbation is related to the density per- 
turbation via the speed of sound. Usually one defines 
the speed of sound in the rest-frame of the fluid where 
the pressure perturbations are assumed to be adiabatic: 
c 2 = (8 P d e/8p de ) rest- The transformation back into the 
Newtonian gauge leads to a correction to 8P de which is 
proportional to the velocity potential, 9 de /k 2 . Here we 
simply parameterise the pressure perturbation as 

SP de _ 2 x ,9P 



Pde 



k 2 



(28) 
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find from Eq. J27|) 




0.1 1 



a 

FIG. 1: The evolution of Q and r\ is presented with varying 
Cp for the case of scale-invariant growth of dark energy per- 
turbations where the anisotropic stress balances the pressure 
support (fc = 1). Here we set Wde = —0.8 as in the sDE 
reference. 



Both, cp and gp can be functions of time and space. 
However, we will consider cp to be constant (52j . 

In the absence of anisotropic stress, a de — 0, dark 
energy perturbations are strongly suppressed on comov- 
ing scales smaller than cp/(aH) because of the negative 
scale-dependent term oc k 2 8Pd e on the rhs of the evolu- 
tion equation (|27|) . Therefore, dark energy perturbations 
become only important on quasi-static scales if cp <C 1. 

The presence of anisotropic stress in the dark energy, 
however, can mitigate this effect by compensating the 
pressure support. The extreme case where Ode = & de 
with 

(1 + Wde)Pde 

leads to a scale-independent growth of dark energy per- 
turbations if gp = 0. In this case the scale-dependent 
suppression term on the rhs of the evolution equation for 
5d e vanishes. More generally we set 

Ode = U<T° de + || , (30) 

such that the constant parameter f a re-introduces scale- 
dependence into the growth of perturbations in case f a ^ 
1 . This effect becomes more apparent when we substitute 
the expressions for the pressure and anisotropic stress 
perturbations into the evolution equation. For the simple 
case when Wde and cp are constants and gp = g a = we 



5 de + (2 + 3C P - 6w de )HSde = 

^(1 + W de )H 2 ^Q m S m + (1 + 3f a Cp)QdeSde 

i + Wde ^ a In a / J 

~(l-U)c P d de (31) 

where we also substituted the source term from Eq. (|26[) . 
For the matter we find 

5 m + 2H8 m = ^H 2 {n m 5 m + (1 + ?>f a c 2 p)Slde5de} ■ (32) 

In the following we will investigate the two cases: scale- 
independent growth of dark energy perturbations under 
the conditions f a = 1, gp = and g a = 0, and scale- 
dependent growth with a deviation from either one of 
those conditions. 

On scales where the quasi-static approximation is ap- 
plicable, the clustering of dark energy is governed by 
the competition between the gravitational infall and the 
pressure support as can be seen from the fc 2 -dependent 
term in the full evolution equation (|27p . However, the 
anisotropic stress counteracts the pressure. For cr de = 
crjL, i.e. f a = 1, the anisotropic stress exactly cancels the 
pressure support. The remaining source term is domi- 
nated by the matter density. Thus, the dark energy clus- 
ters on all scales like dark matter does, and consequently, 
the growth of perturbations becomes scale- independent. 

From the simplified evolution equation (f3i~j) , with con- 
stant Wde and cp, we can read off the influence of the 
sound speed: cp > enhances the Hubble drag, the sec- 
ond term on the lhs of Eq. (f5T|) , and also lowers the source 
term, as the third term in the braces is < — 10 for typical 
values of Wde and all cp > 0. In other words, compared 
to the matter perturbations, dark energy perturbations 
grow weakly compared with dark matter perturbations 
but are not completely suppressed. Therefore we expect 
Q > 1 at late times. Dark matter perturbations do not 
suffer the suppression but are rather enhanced due to the 
dark energy anisotropic stress adding to the source term, 
as seen from Eq. (f52"| . The clustering condition of cp <C 1 
is relaxed compared to the case with no anisotropic stress 
and we find dark energy perturbations cluster on all 
scales even for cp = 1. In Fig. Q] we show the late-time 
evolution of Q and 77 in the case of the scale-independent 
growth. With the help of a modified version of the CAMB 
code [47| we plot the mode k = 0.01 hr 1 Mpc, which is 
representative for the quasi-static regime. It is important 
to note that when computing quantities like Q with the 
help of a full Boltzmann-code like CAMB rather than 
within the quasi-static approximation it is necessary to 
use the correct gauge-invariant definition of Q in terms 
of the comoving density perturbations A, rather than 
the Newtonian gauge quantities <5j. Omitting this leads 
to artificial scale-dependence because it takes time for 
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FIG. 2: The scale-dependence of Q and 77 is shown for different 
values of f a . We include the limiting cases of no anisotropic 
stress [f a = 0, dash-dotted) and scale-independent growth 
of dark energy perturbations = 1, dotted). Here we set 
Wde = —0.8 and c% — 1 as in the sDE reference. 



modes to relax to their late-time asymptotic behaviour 
after entering the horizon. 

Q is enhanced for smaller cp. The anisotropic stress 
that is needed to cancel the pressure support, a<ie = &d el 
leads to r\ < 1 as can be directly seen from Eq. (f!H)) . If 
cp = there is no anisotropic stress, but also there is 
no scale-dependent suppression of the source term in the 
growth equation for dark energy perturbations Eq. ([31]) . 
Therefore dark energy perturbations cluster on all scales 
like dark matter does. Though the growth of dark energy 
perturbations is suppressed compared with dark mat- 
ter perturbations because of a larger Hubble drag and 
a scale-independent suppression of the source term. 

If the anisotropic stress does not exactly cancel the 
scale-dependent pressure support in the dark energy 
growth equation (f!?T)) . i.e. j a ^ 1, then the growth be- 
comes scale-dependent. The term oc k 2 on the rhs of 
the dark energy growth equation becomes important on 
quasi-static scales. For f a > 1 it acts as an instability 
and S c ie would blow up on small scales. For f a < 1 we 
can define the scale 



aH 



(33) 



which acts as a reduced sound horizon such that dark 
energy only clusters on scales larger than X a (as opposed 
to the scale-invariant case where dark energy clusters on 
all scales equally.) A slight departure of f a from unity 
therefore effectively decreases the sound horizon as com- 
pared to the sDE. Therefore we expect dark energy to 
cluster on smaller scales than in the case where there is 



no anisotropic stress. As a result, Q and 77 become scale- 
dependent and do approach unity on scales smaller than 
A CT , see Fig. [5] In addition to this case, scale-dependent 
Q and 77 are generated by gp 7^ or g a ^ 0, as the 
scale-dependent effects from the pressure perturbation 
or anisotropic stress cause the density contrast to evolve 
in a scale-dependent way. We leave the investigation of 
these effects for future work. 



V. INTERACTING DARK ENERGY 

Dark matter is currently only detected via its grav- 
itational effects, and there is an unavoidable degener- 
acy between dark matter and dark energy within GR. 
There could be a hidden non-gravitational coupling be- 
tween dark matter and dark energy. In these interacting 
dark energy (IDE) models the continuity equations for 
cold dark matter and dark energy are modified by the 
coupling C such that 



Pde 



p c + 3H Pc = -C, (34) 
3H(l + w de ) Pde = C. (35) 



where w ( i e = Pde/ Pde- C is the rate of the energy transfer, 
so that C > 0(< 0) implies that energy is transferred 
from dark matter to dark energy (from dark energy to 
dark matter). In the following, we consider two types 
of interactions and identify Q and 77 in these models. 
In this section, we ignore the clustering of dark energy 
studied in the previous section assuming that the sound 
speed is one and there is no anisotropic stress. Also we 
assume baryons are not coupled to dark energy as there 
are strong constraints from laboratory tests. 



A. Model I: C = T c p c 

First we consider the interaction given by C = T c 
covariant form of the interaction can be found in [48H5C 
where it was shown that for this form of interaction there 
are no modifications to the dynamics of quasi-static per- 
turbations. In fact, in the Newtonian regime the only sig- 
nal of the dark sector interaction in structure formation 
to linear order is via the modification of the background 
expansion history. 

The energy conservation and the Euler equation for 
density and velocity perturbations in the dark matter 
and the Poisson equation are given by 



0. 



H6 c + —$ = 0, 



(36) 
(37) 



fc 2 $ = 4irG N a 2 (p c 5 c + Pb 5 b ). (38) 

These equations have precisely the same form as those 
for the non-interacting case. 
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FIG. 3: The time evolution of Q is presented for two different models of IDE. The left-hand side plot corresponds to Model I: 
C — Y c p c with varying T c (1/Mpc). The right-hand side plot corresponds to Model II: C = Tdepde with varying Tdc (1/Mpc). 



However, a non-trivial Q is induced by the non- 
adiabatic scaling of p c such that p m = p c + pf, ^ p s m 
although Pb = p%- The Poisson equation is written as 



47rGAra 2 (p b + p c )8 n 



where we used the fact that S m = 5), = 5 C 
this with the definition of Q in Eq. ([5]) 



fc 2 $ = 4irG N a 2 Q(a, k)p s m 6„ 



we derive 



Q(a, k) 



Pb + Pc 



(39) 
Comparing 

(40) 
(41) 



In Fig. [3] we plot the parameter Q(a,k) from Eq. (14T1) . 
with different values of the interaction rate T r . 



B. Model II: C = Y de p de 

Next, we consider the interaction given by C = YdePde. ■ 
A covariant form of the interaction can be found in j35| . 
The Euler equation for the velocity perturbation is the 
same as in the previous model 



k 2 

H9 C +—<I> 
a 



0, 



(42) 



but the energy conservation equation for the dark matter 
density contrast now acquires a source term due to the 
coupling with the dark energy HH : 



1 de O c 



(43) 



This fact leads to a linear bias between dark matter and 
baryon over-densities. The source term leads to a modi- 
fied growth equation of dark matter: 



4ttG< p b S b 



x<2-3w 



2 Tde Pde 

1 H 

3a H p c 



H 



Pde 
Pc 



PeSc}- (44) 



In this model, there is an ambiguity in the definition 
of Q due to the bias between baryon and dark matter 
over-densities. The baryon over-densities obey the same 
growth equation as that for matter perturbations in the 
sDE reference model 



5 b + 2H5 b - 4irG(p b 5 b + p c S c ) = . 



(45) 



Therefore we define the function Q with respect to the 
baryons 

k 2 ^ = ^G N a 2 Q{a,k)p s m 8b, (46) 
such that the bias is absorbed into Q. Hence we have 

PbS b + p c S c 



Q(a, k) 



Pm S b 



(47) 



The reason why we chose to define Q in terms of 5b is 
that we might be able to use the energy conservation 
equation of the baryons to measure the baryon density 
contrast from peculiar velocities of galaxies, which are 
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likely to follow the baryon peculiar velocities. Dark mat- 
ter over-densities are subject to bias and we would not be 
able to measure them directly. The above expression for 
the induced Q in case II is also applicable to other kinds 
of interaction models, e.g. the one proposed by Amcndola 
[35j where the interaction is described by a quintessence 
field coupled to dark matter. In Fig. [3] we plot the pa- 
rameter Q(a, k) from Eq. (|47l) for different values of the 
interaction rate T^e- 



VI. TRAJECTORIES ON THE £ AND p, PLANE 

In the previous three sections we examined predictions 
for Q and 77 in MG models and GR models with exotic 
dark energy As discussed in section III, MG models char- 
acterised by the BD theory in the quasi-static approxi- 
mation have a distinct path on the plane of (£, (i). In 
this section, we summarise our predictions in this plane. 
For simplicity we only consider the case where E and (i 
are scale independent. 

Let us consider the simplest massless case in the MG 
models described by the BD theory. In this case, as is 
seen from Eq. (fb4"|) . Q and ?/ is expressed in terms of 
a single parameter wbd, i-C- both are not independent 
quantities. If we project this constraint on the E and /1 
plane, we get 



E(fc, a) = 1, [i(k, a) 



2(2 + o;bd) 
3 + 2wbd 



(48) 



and interestingly, E(fc, a) = 1 is the same as GR with 
varying \x only. Note that [i and E are scale independent, 
as the massless case is considered. As is shown in Fig. |H 
the MG models described by the BD theory trace a path 
along the direction of [i at fixed £ = 1. 

In clustering dark energy (cDE) models, it is difficult 
to give physical models for dark energy perturbations. 
We have considered a simple toy model where the pres- 
sure perturbation is modeled as in Eq. ([2"5]) . If we impose 
scale- independence of [l and £, we find g p = and that 
the anisotropic stress should be related to the pressure 
perturbation via Ode = 3Pde/{l + w<ie)pde- Then we are 
left with only one free parameter, either the pressure per- 
turbation 8Pd e or the anisotropic stress Ude, thus again 
there is a unique path on the (£, /1) plane. Of course, 
it is possible to vary SPd e and Ode fully independently, 
but then the scale independence of £ and [i is lost, or 
absolute fine-tuning. 

In interacting dark energy (IDE) models, no anisotrop- 
ic stress is introduced with the assumption that dark en- 
ergy is smooth and, thus, only Q is non-trivial. There- 
fore, the condition 77 = 1 reduces the degrees of freedom 
again to only one such that we have the relation 



£ = fj, . 



(49) 




This turns out to be quite distinct from the trajectory of 
the MG models as is shown in Fig. @] 



S(a=l) 

FIG. 4: Trajectories on E and /i plane of BD type MG models 
(solid curve) , clumping dark energy (dotted curve) and interacting 
dark energy (dash curve). 



These examples suggest that if we can measure £ and 
fi from observations, the path on the (£, [i) plane en- 
ables us to identify the underlying physics of the cosmic 
acceleration. 



VII. CONCLUSIONS 

In this paper, we proposed to paramcterise the rela- 
tion between the lensing potential and the matter over- 
densities, £, and the dynamic relation between the New- 
tonian potential and the matter over-densities, [l, which 
enable us to characterise theoretical models and constrain 
them with observations. If dark energy is described by 
a perfect fluid that is homogeneous (smooth) on sub- 
horizon scales these parameters are trivial, i.e. £ = I 
and [1=1. We showed that E and [i can depart from 
unity in some theoretical models; such as modified grav- 
ity models, interacting dark energy models and clustering 
dark energy models. Interestingly, both parameters are 
related to each other in an unique way depending on the 
underlying theory: 

• With the assumption of the scale-independent evo- 
lution of perturbations, E and [l in Brans-Dickc 
type MG models are described by a single vari- 
able, wbdi which leads to a specific trajectory with 
E = 1. This comes from the fact that there is no 
coupling between photons and the BD scalar field. 

• In clustering dark energy the scale-independence 
of E and [i and the simple assumption that 5Pd e oc 
5pde lead to a constraint equation between the pres- 
sure perturbation and anisotropy stress. Therefore, 
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both £ and (i arc determined by a single variable, 
and consequently there is a unique trajectory in the 
(£, fi) plane. 

• For interacting dark energy models, as there is no 
anisotropic stress induced by interactions, we have 
a simple relation, £ = /i. Nevertheless, a non- 
trivial £ is induced by the non-adiabatic scaling 
of the background dark matter density. 

Those investigations enrich our understanding of large 
scale structure formation, and provide theoretical priors 



on the modified growth parameter space spanned by £ 
and \i. 
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